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ABSTRACT 


This thesis contains discussions of a number of points which arose when 
the author was studying the "paramagnetic resonance line shape problem". 
The so-called moment method is discussed, and a new derivation of the 
moments of the line shape function is given. Single-spin operators are 
introduced which simplify the calculation of these moments. The Green's 
function technique, as applied to this problem, and the decoupling 
approximations associated with the technique, are looked at from the point 
of view of reliability and complexity. As a test of the reliability of any 
decoupling, a theorem concerning the moments of a line shape arising from 
such a decoupling is discussed and proved. The Green's function technique 
is applied to the case of the one-dimensional Ising model with spin %, 
where no decoupling of the hierarchy of Green's function equations is 
necessary. A method of calculating thermal averages for this case, using 


difference equations, is given, 
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CHAPTER I 


Introduction and Summary 


This thesis contains discussions of a number of points which arose when 
the author was studying the "paramagnetic resonance line shape problem." This 
problem is one which has long defied attempts at solution and we do not succeed 
in solving it. However, some points arose which it was thought were, in 
themselves, worthy of discussion. The motivation for tackling the problem 
was that with the recently introduced techniques of double-time temperature- 
dependent Green's functions, some hope of solution might be entertained. 

Let us first briefly review the problem (details are given in Chapter It). 
The physical system consists of a crystalline sample containing N identical 
paramagnetic unit systems or "spins". This system is acted on by a constant 
magnetic field H.. As a result the system has a certain number of eigenstates, 
and transitions between these eigenstates can be induced by a high-frequency 
oscillating magnetic field H(t) = H, cos wt. The resulting absorption of 
the high-frequency radiation is referred to as paramagnetic resonance absorption 
and the problem is to calculate, from first principles, the shape of one of 
the absorption lines. In particular one wishes to know how this shape changes 
with temperature. 

The problem that has been studied is that in which the temperature- 
dependence of the line shape arises solely from the average occupation numbers 
of the various energy levels of the system. In the absence of interaction 
between the unit spin systems, these energy levels would be, except for the 
lowest and the highest energy levels, highly degenerate. The spin-spin 
interactions then remove these degeneracies to a certain degree. There 
exists, therefore, no natural relaxation mechanism that would yield a 


continuous line shape: the line shape function consists of a series of delta- 


functions. One has basically to carry out a quantum-mechanical calculation to 
find the energy levels, and then to follow this with a statistical-mechanical 
treatment to find occupation probabilities. With the inclusion of other 
thermal effects (e.g. lattice vibrations), one would be led to broadened 
temperature-dependent energy levels and perhaps a continuous line shape. We 
do not discuss this aspect of the problem. 

The Hamiltonian of the system consists of three terms: (a) a spin 
Hamiltonian which is the sum of the Hamiltonians of the individual paramagnetic 
units in the presence of the applied constant magnetic field and the 
crystalline field; (b) a spin-spin interaction which consists of a dipole- 
dipole interaction and an exchange interaction; (c) the interaction with the 
applied oscillating magnetic field. With no crystalline field present, the 
energy levels of the individual spins are equally spaced, with spacing hw. 
say. The absorption spectrum then consists of a set of lines at frequencies 


0, Wo» 2W., 3W., ... and to single out the primary line at Wo one does a 


a” *'e 

truncation of the spin-spin interaction (Van Vleck, 1948). If a crystalline 
field is present and the effective spin is greater than one-half, then in 
general the individual spin levels for part (a) of the Hamiltonian will not 
be equally spaced. One then obtains a set of primary lines (with their 
associated secondary ones). A truncation of part (c) of the Hamiltonian is 
then required (Kambe and Usui, 1952) to single out a particular primary line. 
This latter truncation has been performed (Pryce and Stevens, 1950); McMillan 
and Opechowski, 1960) using projection operators which pick out eigenstates 
of the full Hamiltonian of part (a). It is possible to carry out the second 
truncation (Kambe and Usui) using operators defined in the representation of 


single-spin eigenfunctions. These latter operators are introduced in 


Chapter II of this thesis. 


The authors mentioned above (Van Vleck, Kambe and Usui, McMillan and 
Opechowski) did not attempt to calculate the line shape but contented themselves 
with calculating its low moments (first and second). This is referred to as 
the moment method, and little can be said from these moments concerning 
the line shape itself. We show in Chapter IV that our single-spin operators 
simplify the calculation of these moments. 

In recent years, the technique of double-time temperature-dependent 
Green's functions has been used to calculate line shapes directly, and the 
motivation of this work was to find if this technique was of value in the 
paramagnetic resonance line shape problem. As a preliminary exercise, the 
line shape function for the one-dimensional Ising model with spin % and 
nearest-neighbour interactions was studied and is discussed in Chapter III. 

In this case no decoupling of the hierarchy of Green's function equations is 
necessary. 

The technique has been applied by Tanita and Tanaka (1963) to the case 
in which no crystalline field is present. Their decoupling procedure is 
designed to make the resulting closed set of equations as tractable as 
possible. Our attempt at decoupling is designed to be as reliable as possible, 
but unfortunately gives rise to a set of equations which we have not been 
able to solve. As a test of the reliability of any decoupling, a theorem 
concerning the moments of a line shape arising from such a decoupling was 
used, and is discussed and proved in Chapter V. Decoupling approximations are 
also discussed in Chapter V. The result of our decoupling is given in Chapter 
VI and compared with that of Tonita and Tanaka (1963). 

These decouplings of course lead to a set of delta-functions for the 
line shape, and Tanita and Tanaka derived a continuous line shape by 


introducing rather arbitrarily a smearing of delta-functions. This aspect 


of the problem is also discussed in Chapter VI where the general Ising 
model is used as an illustrative example. 

For a treatment of the problem in which the individual spin energy 
levels are not equidistant, the individual spin operators previously mentioned 
were used; it was for this problem that they were originally introduced. This 
problem proved intractable and the reasons for this are discussed, also in 


Chapter VI. 


CHAPTER II 


General Formalism 


1. System Hamiltonian 


Our physical system consists of a crystalline sample containing N 
identical paramagnetic electronic unit systems, or "spins", possessing magnetic 
moments. It is assumed that the sample has been placed in a constant magnetic 
field H. and a high frequency oscillating magnetic field H(t) = Hy cos ft) 
(0) 
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perpendicular to H.. Each spin is defined by a "spin-Hamiltonian" H 


Boe 
which incorporates the effect of the external constant magnetic field and of 
any arbitrary crystalline electric field that may be present. The oscillating 
magnetic field induces transitions between the energy states of the spin system, 
and experimentally, we find that energy is absorbed from the oscillating 
field. The "line shape problem" consists of finding the power absorbed by 
the spin system as an analytic function of ©) , the frequency of the 
oscillating field. 

In addition to the interaction of the spins with the two magnetic fields 
and with the crystalline electrostatic field, there are weak interactions 
within the spin system. Were these latter interactions absent, the line=-shape 
would consist only of widely separated delta-function peaks; these "spin-spin 
interactions" serve to "broaden" the lines, giving them their shape (i.e., 
closely spaced delta-functions of differing strengths occur, and they are 
not experimentally resolvable). 

Mention must also be made of interactions between the spin system and 
the lattice. We take it that the temperature is low enough so that the 
effects of lattice vibrations can always be neglected. 


The Hamiltonian describing our paramagnetic system may be written 


Woosh HO (2-1) 
where 
eH = WHO oa (2=2) 
Ho > WM) (2-3) 
HW > HO | (2-4) 
WY = ZH is the Zeeman, or unperturbed 


Hamiltonian with i labelling the (2-5) 
lattice sites 


Hye Wis the spin-spin interaction Hamiltonian 


> 
d including dipole-dipole and isotropic exchange (2-6) 


interactions 


) 
ph - M-Hicoast s the external perturbation, with M the magnetic 
moment of the spin system. 


(2-3) and (2-4) are weak-interaction assumptions. 

We assume that na possesses axial symmetry, the axis being the direction 
of He, which we choose to be the z-axis; and that the eigenstates of HW are 
also eigenstates of the operator corresponding to the z-component of the spin. 

Before discussing the line-shape function 7 (w) or its moments, we 
must modify the Hamiltonian in two ways. The term oe in the Hamiltonian 
must be "truncated" in order to eliminate the weak "secondary" lines from 
The argument for this has been given by Van Vleck (1948). This truncation 
consists in keeping only that part of av (we will call it tT ) which 


commutes with W ; i.e. 


[pH] = 0 (2-8) 


We further define 


A ” mph > vs WI) =Wy = aka ue (2-9) 


Then [a WW) = 0 (2-10) 


Also, when a crystalline electric field is present, and, thus, the energy 
levels of an unperturbed spin are not in general equidistant, Wo must be 
modified in such a way as to cause transitions between only two levels of a 
spin. This must be done in order to ensure that other "primary" lines of 
comparable intensity to the one under consideration, are not included in the 
calculations. The necessity for this further truncation was pointed out by 
Kambe and Usui (1952). We let Zo represent the modified #”, and we 
put a circumflex over its associated operators. 

When no crystalline field is present, and He is in the (-Z)-direction, 


—_— 


the Hamiltonian Wma. is given by equations (2-9), (2-5), and (2-6) 


with 
re 
) be (2-11) 
Q\ >= 
Ce epe Ay SoG + By SPS (2-12) 
or 
wa —* sat a =p = 
(We), af, > Ai [Ss S +S) S ) + Cy af JS; (2-12b) 
where 
“iP aN ay 
_ - ~ a a = 2 a 
Ay a Ai + J / fey (34;~1), (2-13) 


a he (agte at) say (34-1) 


vj is the direction cosine of the line of length ri joining spin i to 


spin j relative to the z-axis; A,. are the exchange interaction constants; 


t 


g,. and g, are the g-factors perpendicular and parallel to the symmetry axis; 


- is the Bohr magneton; 


and (2-14) 


W). = 4 lH. > 0 
is the unperturbed absorption frequency. 

Where a crystalline electric field is present, the unperturbed sping have, 
in general, non-equidistant energy levels. We consider here only those special 
single-spin Hamiltonians HD which have the property that, if IX? are 
the simultaneous eigenstates of a and af? with eigenvalues A and -1y 


respectively, i.e. 


6*IND> = dAIAZ, 


‘ . (2-15) 


WP 1X2, = - TG IAD, (2-16) 


then 


T), 


+| ~ 1) 9 0 il ad 
Equation (2-17) signifies that in the extension from the no-crystalline field 
case to the crystalline field case, the order of the unperturbed energy levels, 
which in both cases are labelled by the eigenvalues of $.°, remains unchanged. 


The statement about non-equidistant energy levels than means that it is not 


generally true that 


“Pe Oy * fom 4g (2-18) 


when d “ x , 


The “Wi. given by equation (2-12) does not any longer commute with 
(0) >) 
eH” = 2 4 where "> is defined by (2-16). It must be truncated 


still further. The truncation can easily be done if one defines operators 


Ay by the relation 
<x IAS IXZ a by d ‘ S\~ 2 aaa 


From definition (2-19), equation (2-16), and the well-known commutation 
+ 
relations satisfied by the ordinary ladder operators oF, » the following 


equalities can be derived and are found to be useful: 


i 
[45 $f Jab * F 1 A (2-20) 


LS 4°) - ‘ar he (2-21) 


Dy Ay * Dy ae Aix Aran (2-22a) 
+ 29 : +. 
Pex Ay d* = }: a! Ax» Ay 01 (2-22b) 
+ = +, 
Ai) Diy = Dy el a: > Aj)+! (2-22c) 
As) Ass = by “del Axx fas ee (2-22d) 
With the use of equations (2-20) to (2-22), we find for the crystalline field 
case ] - at :) 
iv) +A. S_ A Ay + A; 
a ty DX 41 A bX? 
(i) Yomeg ‘ 


-2 
a C; irs” (2-23) 
where ware means summation over all ordered pairs (A" 1’) 


# I am indebted to Dr. K. Nishikawa for the suggested use of these operators. 
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satisfying 


yn - Ty, = Vy ee oe? ; (2-24) 


Using 
(AS Y* = Appa (2-25) 


which follows from definition (2-19), one can verify that Hie is 
Hermitian. (The asterisk will always be used to denote the Hermitian 
conjugate). 

A Hamiltonian which we consider at some length in a later chapter is 
that of the Ising model. This Hamiltonian was first proposed by Ising (1925) 
to describe certain ferromagnetic systems. It turned out not to be valid 
from a physical point of view, in that it does not describe, with any 
accuracy, the interactions taking place within any known spin system. However, 
it is invaluable from a mathematical point of view, in that it is a 
specialization, of equation (2-12), which is easy to handle. 

The model Hamiltonian is given by (2-1), (2-2), (2-5) and (2-6) with 


_— 


un am Wy S; (2-26) 


— Zpr 
Ca a é, S; S (2-27) 


where the ey are the exchange interaction constants. It is seen that 
(2-27) is just the general paramagnetic resonance Hamiltonian (2-23) with the 


a i and the Ci formally put equal to zero and ~ey respectively. 


In many instances a further specialization is useful, namely that of 
putting 


é 4 nearest neighbours 
4. (2-28) 


‘J fa) otherwise 
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The Hamiltonian with this restriction on the Ey will be called simply the 
"Ising model Hamiltonian". When this restriction is removed, we will 
specifically refer to the "general" Ising model Hamiltonian. 

These are the only Hamiltonians that will be considered in this thesis. 
It remains now to exhibit 1, for these various cases. 


We have, from equation (2-7), 


WO = - MH, (t) ; (2-29) 


as the oscillating magnetic field is taken to be in the x-direction. 

No truncation of Mx is needed for the case of equidistant unperturbed 
energy levels, and so in the absence of a crystalline electric field (and this 
includes the Ising model cases), Mm, = My . However, for non-equidistant 


energy levels, i.e. in the presence of a crystalline field, one must use 


(ha). = ei. 2, a. GN + Aine ) (2-30) 
= ot 


where — means summation over all X such that 


haa Ty ew, «OAS &,) . (2-31) 


and W, is the unperturbed frequency of the primary line under consideration, 
not necessarily given by (2-14), but always positive. 

It is of practical interest to know that H),, ef, is a special case 
of He, » and that (Md) re ef, iS a special case of (Mde, . To show 


this, one needs the help of the relations 


a9 a 
ZA) = S; (2-32) 


t 


A&G, th, (for no of.) (2-33) 
A ye Gr, obh eM i -f) (2-34) 


where in (2-34) it is assumed, naturally, that ae and ee exist (see 
equation (2-24)). (2-34) is a consequence of (2-33) and the definitions of the 
sets @ and G, . Equation (2-32) follows immediately from the definition 
(2-19) of the 4,5 

In developing our formalism for the line shape function, then, only the 
Hamiltonian describing the crystalline field case need be taken into 


consideration. 


2. General Formalism of Absorption Line Shape. 


We start from the well-known relation between the absorption constant 


o~ (>) and the magnetic susceptibility X (wd): 
ride eont Win X(v) (Wo), (2-35) 


In the Kubo-type formulation (Kubo, 1957; O'Rourke, 1957)  (w) is written 
as 


X (w) = 4 gin (ae e ei <IMx (t) M.J> (2-36) 


where we have used the following notation: 


My : the x-component of the total magnetization operator 


oie 62° Oar Fhe (2-37) 
‘4+ = will henceforth be called the Hamiltonian of 


the system 


5 __ , -K/k; 
Cal a O'Ta gat Q= 5 (e if) oe 


) 
k, Boltzmann's constant, 
T =: temperature of the system at the time t = - od. 
The following assumptions were involved in these formulas: 


(1) The external magnetic field H(t) has been applied in the x-direction. 


(2) The response is looked at in. the x-direction. 

(3) Hy (t) has been applied adiabatically starting from the distant past: - 
lim Hy (t) = 0. 

(4) The system was in thermal equilibrium at t = - co. 


(5) The Hamiltonian associated with the oscillating field is of the form 


Po = - MH @) (2-39) 


(6) The oscillating field is sufficiently weak that its effect can be 
treated in first order perturbation theory. 
One normally defines the paramagnetic resonance line shape function to 


be not 7 (2) » as given by (2-35), but g tu) as given by 


4 (w= const. Lm X(W) (W>0) (2-35a) 


where the constant is such that 4(\)) is normalized to unity. (Cf.McMillan 
and Opechowski (1960)). 

It is of advantage to keep only that part of X (w) as defined by (2-36) 
which vanishes for w<0 . Employing equations (2-8), (2-9), (2-37), (2-30), 
(2-21), and (2-31), and the well-known identity 


P 
a a? se = Lats (2-40) 
where fF, a [... [F, 6),6), = 6] { F, G}'= Y. (2-41) 
one has P Gs 


Ms, ¢ (-t) = eri HEY e Ree at (He6 peye AN etfs Rut 
= einer ( en HOt eb WHE) o THE 


oni #" i f) iyo 
= (Agu € “(eee ’ttpamave' "Je t 


et swe ° -thy, - YL a - 
“Ads & eal (&= fel ar, 5, ep (es4e) 
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With the use of the identity (2-43) which follows from the cyclic properties 


of the trace for any operators A and B 
<Casey>= <CACd,WD> (2-43) 


and the relation 


i a(S +S) (2-44) 
+ Age 
equation (2-36) can be rewritten as follows: 


a jie stlewd t-et (MS ELE, Gre git, . ae 


Pian SE poen We, b> 
e i (a \t- et ae a “e- Ye 
+h ey ite ea ve ae Palos 
*  _i(oent-e0 gthit Epi 
+i Gin (ac é€ (444) aE ste dyn & i, 
GFE? o* ° vies ¢ wt wep 
A : - = 
+1 Yn (de en ilotheje-eL a f iS Ai @ | 
Fé o* : fpe= i ‘eG (2-45) 
First, notice that in the second integral, See 4 may be 
replaced by Ze dia (see Appendix I). 


As the system is definitely quantized along the z-axis, the first and 
fourth integrals vanish automatically. Moreover, the third integral provides 
a negligible contribution to X(wW) for -W > 0; this follows from the weak- 
interaction assumption:(2-3) if one puts }) =O, the integral gives a 
delta-function peak at W=- We  . Inclusion of >” — serves to 
introduce many other delta-function paaks into the region around -We. , with 
very little overlap into the region of positive ™ . By a similar argument, 
the second integral is seen to provide a negligible contribution to X Ww) 


the region of negative W .# 


# This argument is equivalent to that of McMillan and Opechowski Seti They 
argue that if W” is not too large, Ex< E, implies Exai < k except ing 
a negligible number of pairs of states Ix:),/4,R) where p¥) |4,k¢ See ae 
and 4 /#,k)= Eph |A,k) (ke 1,294 . 29,» the degeneracy of E. 
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Thus, if the magnetic susceptibility is redefined in terms of the second 


integral alone, i.e. 


Xl) = Ce.) eo * Lain fire Pal & Coke A1y4 Ct, a = 4:4 iy, 


then (2-35) may be rewritten without the condition" w>9o ". Further, one 


— has 


fi D(a) X Ce) a -[ 0) X (ao) old ae 


where ) (1) is any function of W 

The function X (w) as given by (2-46) is closely related to the 
temperature-dependent double-time Green's function. .These Green's functions 
were first introduced by Bogoliubov et al. (Bogoliubov and Tyablikov, 1959; 
Zubarev, 1960). The Green's function €A| BY, of the operators A and B is 


a two-branch analytic function of E defined everywhere outside the real axis 


A By | om a. - 
K¢ E. 7 ~Lpr pty? cama, f —— a 48) 


By comparison of (2-46) with (2-48) one can immediately see the relation 
& - + 
em SE 4h, (Sa Pew! 
x (0) ar (gu) Lom, < Gt mre th Manele aah 
a L} i eA) > at 
= = 2m 89.) <== iy 24, ll ila 
where we have used only the identity a 
Note that the above Green's function is of the form CAIN (by (2-25)). 
One can show that the Green's function ¢ AlBDe as defined by (2-48) 


is analytic both in the upper and lower half-planes and that it is found from 


the hierarchy of wath 


e Kaisz = ae <laer <[A,¥] Is 
E SAMIR = <CL4,#2,ap + €(LA, 4], HIB 


2-50) 
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Our problem is then to solve this hierarchy with the condition that the 
solution be analytic in-the domain: for which the Green's function is. defined, 
that is, outside the real axis. .The.actual procedure of solving this equation 
approximately will be discussed in Chapter V. For further details of the 


Green's function method, we refer to Zubarev's article. 


Lf 


CHAPTER III 


Green's Functions..and.the. Isin Mode1* 


In some simple cases it is found that the hierarchy of Green's Function 
equations (2-50) decouples. automatically after a finite number of steps. One 
such case is the problem of finding the amplitudes and frequencies of the three 
component lines in the absorption band of the one-dimensional Ising model, 
consisting of an arbitrary number N of identical spins arranged in a ring. 

It is assumed that each spin has spin quantum number 43, and that only nearest- 
neighbour interactions occur. 


The Hamiltonian of this system is given by equations (2-26), (2-27), and 


(2-28), i.e. 
= N 2 + 2 = 
aN = - h 2->e € Jesiai Sig Sis (3-1) 
where 
h = We as given by equation (2-14) 


m 
" 


exchange interaction constant 
£ and ke label the lattice sites, with the (N+1)S* site identical with 
the Ist. 

The line shape function g (w) is again defined by equation (2-35a), 
which, through equation (2-49) specialized to the no-crystalline field case, 


leads to 


4 (13) al. le Ln Gs hgske (3-2) 


# I am indebted to Professor W. Opechowski for his suggestion that this 
problem be considered. 


where 


Cele = = o Clas IS7R (3-3) 


The first three Green's function equations are then 
(E -h)KSEIS “2s AOD Uy +E CSS (SE Ss Se 2, (3-4) 
Eh) KSSH = ENC 6 CG SE+GEY SQ 


(3-5) 


E-h)KSE GGG N= Ge OP Fete God sr gay h, 
FERS (E+) ISR 
= 4236 (<2? + 4<Gtstse>) 
+E KSEE GF) ISSEY 


(3-6) 


The last term in equation (3-6) arises since 


(8 Ghy= (GF agrt 


(3-7) 


is a consequence of the relation 


(4) /4 ($= ifn) . (3-8) 


) 
It is now seen that equations (3-4), (3-5), and (3-6) form a closed set of 


simultaneous equations which can easily be solved for Ge le ; that is, there 
is no need to write any higher-order equations. Since this state of affairs 
appeared without any decoupling becoming necessary, we say that the hierarchy 


decouples automatically at the third stage. The solution is 


KSt1SpR < oe re (<7 - 4<3) tre (<1? - 427 +37) 
(3-9) 


+ pee tiarket 4<37) 


where we have introduced the notation 

e= E-h 

J>= <5) (3-10) 

<a) = Eh ey 

cape S55 SSP 
noting that 43 , <2> » and <3> are independent of i because of 
translational invariance. 


The line shape is then given, with the use of (3-2) and (3-3), by 


Cw) Lin com G, lgew-te 
norm, Es>or norm, 


= = [ (<i = 437) lob) 


a<iy 


+ (460- 2<a4.2<39) 0 Cu- hee) 
(3-11) 


+ (47+ adapt 3.<37) P(w-h- 6) 


which is a sum of three delta-functions at specified frequencies, with 
amplitudes given by. linear combinations of thermal averages. 


It now remains to calculate these thermal averages. Since several 
authors (e.g. Huang, 1963) have derived explicit representations for the 


partition function Z, we will take it that Z is known, and use Huang's 


result 


Zz = (r.)" 7 (\.)" (3-12) 
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where 


has fer [eesh (4h) = [essh* Chih) 2.67% ainh 490) (3-12) 


We immediately get <1? and <2> by direct differentiation from the 


formulas > 
| eee 
ee pM (3-14) 
| re 
7 = NA Me) (3-15) 


The thermal average ¢39 is obtained as a. linear combination of <1) and 


<2) by the use of the following general identity (see, e.g., Zubarev, 1960) 


7) 
-iwlt-t’ 
<BAe’p * ; fave . fea lO te ~ 410.6} ee 
“00 am v | 


Putting t = t’ =0, A= s,’, and.B = S,_ 


F ; in equation (3-16), and using 


the result (3-9), we get 


eo de = ee 5. & (<9? - 1a + 4 <37) 
ef" = | eath-« ~* 


(3-17) 
Eh (6 + 4 day + 437) 
eA th+é) _ | 
or, 
~o= D+D<7 + DL? (3-18) 


where Do» Dis and Do are simple functions of h and € . In the derivation 


of (3-17), use was made of the relation 


ee - Si (g= 4a) (3-19) 


) 
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Equations (3-14), (3-15), and (3-18), when combined with equations 
(3-12), (3-13), and (3-17), give us, when substituted into equation (3-11), 
an explicit representation of the absorption band of the one-dimensional 
Ising model. The result corresponds exactly with that found by Y.Y. Lee 
(1961), who had quite complicated summations to perform. The latter arose 
from his combinatorial approach to the problem, and necessitated the finding 
of suitable generating functions. He simplified the combinatorial problem 
by specializing to the case in which N = 4P -1 (P=1,2,...) and then 
considered his results valid for large N. In our approach there is no 
need for this specialization, as the complicated combinatorial problem is 
avoided. The advantage of our method lies in the identification of the 
amplitudes of the resonance lines with easily recognizable thermal averages, 
and in the ease with which these averages can be calculated, from a knowledge 
of the partition function, through differentiation and further application 
of the Green's function technique. 

The thermal average <37 cannot be obtained directly from Z; nor can 
higher-order thermal averages. A direct method of finding thermal averages, 


using difference equations, is given in Appendix II. 
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CHAPTER IV 


The Moment Method 


1. General Discussion 


In the problem of the general paramagnetic resonance line shape, 
unlike that of the one-dimensional Ising Model, the hierarchy of Green's 
function equations (2-50) does not decouple automatically at an early stage. 
The solution of this problem by the Green's function technique then 
involves looking for a decoupling procedure which gives a closed set of 
equations satisfying the following requirements: it must be (a) simple 
enough to be tractable, (b) complicated enough to yield far more than the 
small number of isolated delta-function peaks obtained for the Ising model, 
and (c) "faithful" enough so that it furnishes a line shape function that is 


resonably close, in some approximation at least, to the experimental ly- 
observed line shape. 


The difficulties inherent in the search for such a suitable decoupling 
procedure are representative of those encountered in other approaches to 
the line shape problem. Consequently, a method has been evolved which 


enables one to discuss the shape of the resonance lines without actually 
having to find an analytic expression for the line shape function. Known 


as the "moment method", this procedure consists in defining a line shape 
function in a formal way, and in calculating its first few moments in 
various approximations. It was introduced by Van Vleck (1948), and used 
by Pryce and Stevens (1950), Kambe and Usui (1952) and McMillan and 
Opechowski (1960). Kambe and Usui derived quite general exact expressions 
for the moments, and these were applied by McMillan and Opechowski to the 


case in which a crystalline electrostatic field is present. In this 


cs 


section, a derivation, different from that of Kambe and Usui, of the 
expressions for the moments is given. 
We shall show that if a normalized line shape function Ina Cw) . 


defined by 


dan (i) = em Kaul) /{ Tm Yonlwrale 


(4-1) 
with ” 
: .! -c(w-lé ag 4-2 
Xeaalw) 4 Figen fa lias <la,2e)]7 
satisfies the conditions that 
A* = 8 (4-3) 
and that 
Deva (w) = Oo ,Ww<sO, (4-4) 


then its pth moment about the origin, Mp » iS given by 


p= <[ fA, ap ay? /<La, Ry pepe: (4-5) 


_ en 


With the use of (2-40), (2-43), (4-3), and ie identity 


<Ca (t),5)) = Le a*t)] > (4-6) 


we can immediately write, from (4-2), 
(.3) 


on Pan YS (i/at) Kos (io) — Wes 
a ae + {at eee ee Fact) B)> 


E>o* 
7 Li 1k (dee or iwt- -clt| 3 <a) | 


(4or) 
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We now make use of the characteristic function M(t), defined by 


y 6g tut 
M(t) = [4 e : 7 BA (w) (4-8) 
-DX& 


which possesses the inverse Fourier transform 


dw) = an [at € 


Repeated differentiation of (4-8) with respect to t, i.e. 


= Lie 


M (+) (4-9) 


wd : 
p ) a. oF ( =< rt 


gives, simply, the Taylor coefficients in the expansion of M(t) about 
t=0, i.e. 


=— Get)? 
M(t) = 2. a (4-11) 


Substituting from (4-11) into (4-9), we get 


a wt Zab P 
Das (yw) ar dte ‘ 2, cory (4-12) 


~@ 
Comparison of (4-12) with the alternate expression for 95 a (>) given by 


(4-1) and (4-7) yields 


ar: = (constant, independent of p) X af Can} a} (4-13) 


The constant is determined by the normalization requirement on Gen (uw) . 
resulting in the required expression, (4-5), for the moments about the 


origin. 


If, further, my can be expressed according to (2-9) as 
W = yO 4 EM 
where (2-10) holds, and also 


[a 4) = WA (4-14) 


then 


_ <Cfa, aH BP 
M. < fo-w'g, taal Pag scab (4-15) 


are the moments about the unperturbed frequency 1 (We call the M 


"central moments". ) 


For 


Pee foPdaa(v)aw = [stn Co) ded 
(es "Cr (w- ar ys: “Gan! (v2) eli) 


i) 


ig 
wc p-+ 
4 2 ‘fie W, M.- 
: (4-16) 


Also, by induction it can be shown, using (4-14) and (2-10), that 


x A} Ss 'c, we L ABP (417) 


Combining with (4-5), we . 
ae on WPM, s Fr Le WwP4CEA, opi B]? 
< LA, BI? 


(4-18) 
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Taking p=0,1,2,... successively, we see that formula (4-15) follows. 


2. Application to Paramagnetic Resonance Line Shape Problem. 


For the paramagnetic resonance line shape problem, one can immediately 


proceed to calculate moments in various approximations by substituting 


As 2 ah a 
and 
B= A = “a Zs, Ajj A‘ (4-20) 


into the expressions (4-15) for the central moments. The expressions (4-19) 
and (4-20) for A and B come from equation (2-46) for xX (w), One takes 
commutators with aaa) (equation (2-23)) according to the rules (2-20) 
and (2-22), and takes traces in the most convenient representation. This 
procedure is straightforward. 

McMillan and Opechowski (1960), starting from formulas first derived 
by Kambe and Usui, were able to cast expressions for the first and second 


central moments into a form suitable for computation. They then applied 
their results to the specialized problem that we have been considering. 


Their only initial assumptions are given by our inequalities (2-3) and (2-4). 
The procedure for their calculation of the first central moment, M)> will 


here be given in outline. They start with 
Ma a Sa (b* CA CH A) non 
| => -~ - 4-2] 
rr (b" [ Me, 4) 
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DL 


= 2xp (- 1 /kT) 


2 + is the component of the magnetic moment of the spin 
system in the direction of the oscillating field, truncated so as to 
single out a particular primary line. 


A 
M, and M_ are defined by 


a : (Aiilm lah) if 1E,-6, |= and By Es bk 
(4,1, 18, k) = { 0 otherwise (4-22) 


(A LIM | k) if le,7E,|= We anol Erk read 


bet,i [M.la,k\= } O otherwise (4-23) 


where | 4 :) satisfies the eigenvalue equations 


Pte i) =e, IJ, c) (i= Party 2 the degeneracy of yd 


AL [Av) = Bet [A i) (4-25) 


They then rewrite equation (4-21) in terms of # and M rather than # , 


a A 


A 
M, M, and M_.  . This is done with the use of projection operators Py 


defined such that 
Pr lak) = Oi IA, k) (4-26) 
With the use of the observation that if + is not too large, then 
Ey < En implies Eu; "4 Ea, k except for a negligible 
number of pairs of states 14 i) ; IA, k) , the definitions (4-22) and 
A “n 
(4-23) of ML and M_ are changed to 


CAL [Mla,k) if £,4-E, = Wd, 
O ‘otherwise (4-27) 


cn 1M. 1,k) = { 


Shei. . 


(4,(19-1A,k) = io sil ia it Rasen “tae (4-28) 


Then ite p 
> = po (4-29) 
Me 5 fe ME (4-30) 
M_ = —_ P, Mh (4-31) 


En -£,= WwW, . (4-32) 
Equation (4-21) becomes 
te 2 2/6 (4-33) 
where 
RHP fs HE 
oe Soe fb “amp, Mb" "a MPM} (4-34) 
pawre (1) 
ri aft (MP, WB a AM pM) (4-35) 


NACo We MEM Boe HRM), 


McMillan and Opechowski make use of an approximation technique used 
by Pryce and Stevens, and also by Kambe and Usui. This technique depends 
on the weak-interaction assumption (2-3). One expands all exponentials 


involving Bs WH. H as follows 


28 
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CMC. S @) ' ' 
eC = ¢. oe ee [i+ Ee mM"C f. ai US, HC,)* tu] 
provided Pon, . hatte =6 , and keeps only the 


lowest-order terms in the square brackets. Here Cy and Cy are any operators 
such that the above provision is satisfied. We call the approximation 
which consists in keeping only the first £ terms in the brackets, the 
ae th approximation of McMillan and Opechowski", and denote it by (MO), 
Evaluation of formula (4-33) in (MO), consists then in replacing 

b * FH Pa by L&« Pa » The analysis is done by McMillan and 
Opechowski for the case when the unperturbed spin has R energy values 
A> Aas ag, «++ ap (the equivalent of our T, ) which are all non- 


degenerate. In this connection, a set Gs is defined such that if 


YT is in G, , there exists an integer +’ such that 


Ane - Aw = We (4-36) 


In the first approximation the term (3 (4-34) can then easily be 


reduced to the form 


(- bre) =. a a [Lil Ma,k) | (4-37) 


In the above form, the calculation of B becomes a combinatorial problem, 


and involves doing summations like 


as L (1,4, + NG yt 1.+ Na Ge) 


Nyt tit Nes M 
en $N ; 
(o ¢v ) . z. | Cr'lM I7).| (v-1) | 
Teg, ? Es 


th 


(4-38) 


where I~), are the unperturbed eigenfunctions of the i~ spin, with 


eigenvalueS G. so that 
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We lr), = Gx |v); 6a 1, Ba PY (4-39) 
ee 

The cause of the combinatorial problem is in the degeneracy of the ca 

The degree of complexity increases when expressions like .c(1) and those 
entering the higher moments are considered, and when approximations higher 
than the first are used. 

The combinatorial problem is avoided for the special case that we 
consider by the introduction of the operators ae defined by (2-19). 
To see how these operators simplify the calculations for our special case, we 
show how the more general case of McMillan and Opechowski is simplified by 
the use of a different kind of projection operator; these projection operators, 
first introduced by Kambe and Usui (1952), are then related, for the special 
case, to the operators an 

We define projection operators Py which project onto eigenstates of hl 
rather than the operators P, which project onto eigenstates of the full 


Zeeman Hamiltonian, by 


p* aed, 3 ie LOF (4-40) 
The relations ce 
ame ‘3 = any’ r. ii 
and 
[A 47 ‘ee (4-42) 


follow from the definition (4-40). 


3] 


* * z 
Note that in these equations A no longer labels eigenvalues of JS) r 
(0) 
but only eigenvalues. of HH: according to equation (4-39). 
A A 
One can then show that the following expressions for M, and M_ satisfy 


the definitions (4-27) and (4-28): 


ware 3 ; (4-43) 
A p - 
M, 22. 2." % iy ea (4-44) 


NN 
where i ee | (4-45) 
One also has 


—, a, 4 J a d + 
ia = Ps an = fh, i Hy be A (4-46) 


where - - means summation over all A, ; dy such that 
A3A4 
3 


ay, a ay, = A, + qQ ) L (4-47) 


A 
One can prove the above statements by considering matrix elements of M,» 


A iin 
M_, and 4" (as given above) between the orthonormal product eigenstates 


J\Ay Ag «. 4) = | d.), x Id, x As), XX [Aw)n (4-48) 


' : ) / , 
Here Id,), » for example, is an eigenstate of me with eigenvalue 
q) - The product eigenstates (4-48) are themselves eigenstates of 
‘ 
) 
(9) = ‘ j 7 = 
ph? 3 ona with eigenvalues ay Pay Foss a) E » and are 


linear combinations of the JA, i) defined by (4-24) and (4-25), for fixed m . 


Let us now evaluate 


B= ri Nall [ Me, 47) ‘ (4-49) 


ge 


the first approximation part of the denominator of M, given by equation (4-33), 


with the use of the Py . Using (4-43) and (4-44), 
i 
Mg (Lb % w7 We “4 / d 
a * a 22 2 2 (FM: 1 BE js Tt, 50) 


From (4-41) and translational ar sk 
B= Wa bem (2,'M: 2. My Pe 
_ pt M: BoM Pes)? 


Let us consider only the first term of (4-51), and expand the trace in the 


(4-51) 


orthonormal product eigenstates (4-48). 


B (i teem) = N Be f Aedve Al TE 
«(P'M: PEM: P!)PAs, Aas, Aw) f 


©) 
[%e 


NS FS TOS ol ene ae) 
psi i) 
/ Qu vi 
% eae: [ Gr} mbr’), | 
(4-52) 
where 
i ae ln ass 


One derives B (second term) in a similar fashion to get 


B=NZ2" (. ope i" ede (4-54) 
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where 


(-r’|m|~) = (r’lM; [7), independent of L (4-55) 


L 


and use has been made of equation (4-36). 

There was no need to solve a combinatorial problem, and result (4-54) 
agrees with McMillan (1959). All other traces can be done in the same manner, 
even in higher approximations. 

For the special case when (a) He has axial symmetry; (b) the external 
magnetic field Hi is directed along the axis of symmetry (Z-axis); (c) the 
order of unperturbed energy levels is preserved in the transition to the 
no-crystalline field case (equation (2-17)); and (d) the relation 


[5 " HE] = © holds, we may write 


My, ie 1 Cand a. Mie, (4-56) 
te AE G, 
A + 
M. * vag oo oe A;) (4-57) 
~ te AEG, 
Baer nT tis . 
= We as given by (2-23) 
ity ( e+. 
+ = 
where A ‘ A; 7 Any » and G, are as defined in Chapter II. A 


direct comparison between equations (4-56), (4-57) and equations (4-43), 
(4-44) cannot be made because of the differences in the definitions of G, and 
9 and in the labelling procedures. However, it may easily be verified that 


G 
Ms as given by equation (4-56), conforms to the definition (4-27), and that M, 
as given by equation (4-57), conforms to definition (4-28). 

It has been checked, for the specialized problem wherein conditions (a) to 
(d) (above) hold, that formula (4-15), in the case of an abborption line from 
the ground state to the first excited state of the unperturbed spin, yields the 


same final results for the first and second central moments. as were obtained by 


McMillan and Opechowski (1960); the calculations were done in the first two 


approximations of McMillan and Opechowski, for the cases of effective spin 


S = 34 and S = 1. 
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CHAPTER V 


The Moment Theorem 


1. Relation of Moments to Terms in Hierarchy of Equations. 


It is interesting to take note of the identification between the raw 


moments (moments about the origin) and the first terms on the R.H.S. of the 
Green's function hierarchy of equations (2-50). The latter may be written 


4 i 

E <fa my G7. = s= <Ciaay apr <ha Wt lax (5-1) 
p = GF fra» 

Comparing with the expression (4-5) for /+p » one sees that, apart from a 


factor ar /<CA, 8] the identification is complete. It is assumed that 
conditions (4-3) and (4-4) are fulfilled. 
Moreover, if condition (4-14) holds, the hierarchy of equations (5-1) 


is easily cast into the form 
E-W)< 6A, WP B= se <(LAWV EP AH} BZ (5.0) 
oe Ooleev.<eee 


It is in this form that the Green's function hierarchy of equations will be 
used in the remainder of this thesis. Apart from the same factor aT Kia, 
one identifies the first terms on the R.H.S. of the latter equations with the 
central moments, Mp (equation (4-15)). 

One may ask whether this identification can be made directly from 


equation (5-2) and the relation 
Lin Tm €AIBPe = o-te 
Fa 10) = 
sab [old In CAIN -w-re 


(5-3) 
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between the Green's function €AIBE and the line shape function 9,.¥) A 
proof that this can be done is given in Appendix III. 

This identification is, of course, valueless if one wants only to calculate 
exact moments; the Green's function technique need not be considered at all 
for this purpose. However, very often the Green's function technique, coupled 
with a certain decoupling procedure, is used to derive an approximate analytic 
expression for the line shape function, and it is asked what the associated 
moments are and how closely they approximate the true moments. The integrations 
involved are often extremely difficult to carry out. A theorem, closely 
related to the aforementioned identification, can be proved which provides a 
relatively simple method for finding the approximate moments, and eliminates 
the necessity for performing any integrations. Before this theorem is stated 
and proved, however, an outline of decoupling procedures as they are normally 


followed, will be given. 


2. Decoupling Procedures 


First, rewrite (5-2) as follows: 


Bg PFT + Gia (5-2a) 
where we have replaced S"*Ws ty “CG (5-4a) 
Cla Hei By by G, (5-4b) 
ma P on am 
and — <({An 3) by | e (5-4c) 


A method one might follow would be to write the first q equations as they 
stand, and in the qtlst to replace Goo by a linear combination of the gtl 
preceding Green's functions, the coefficients (usually thermal averages) 


being independent of E; i.e. 


3/ 


€.G, = To + G, (5-5) 
& G, = Fa . Oy (5-5a,) 

ee a 7 (5-5a_) 
e. Or, g-. + Or %q 


eC. oF = Ta a a Bede Gy (5-5aq., ) 


We call this "decoupling at the (qt+1l)st stage". A closed set of linear 
equations is obtained which must be solved for G,, » which is simply related 
to Gan (1) . (It will be obvious throughout when G, stands for the exact 
re and when it is used to mean the approximate value obtained on solution 
of the decoupled set of equations.) 
This method is fourid to be unsatisfactory in most cases. Usually G is 


A) 
made up of more elementary Green's functions Ge such that 


= €22 45/2 = Aji d (5-72) 


and 


aA, chk \’ dat ai Ne | (5-76) 


Let us write a separate hierarchy of equations for each of the elementary 


Green's functions Ot “ 
eh) my?) A) 
e a e he eS" 


(5-8a,) 
(A) wegietl) (A) 
Agee Ty  & (5-8a,) 
Fj | t ! 
(%) | ) 
e Gs Po vs + Go (5- -84q41) 
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Equations (5-8) represent n hierarchies of equations, each hierarchy labelled 
by a different A . They differ from equations (5-2a) in that different A 
A) 

and B have been used in the defining equation (5-2). The T_ are the 

: ; : Cy) (A) ) 
thermal averages associated with the Green's functions oe The G, (r+9 
are not meant to be elementary Green's functions, but are defined by equations 

(A) 
(5-8), through the given G, and equations (5-2). (There is no difficulty 
.) , 

in picking out which terms are to be grouped with the a and which with 


(a) : a ae 7 
the Go. in writing equations (5-8); all the G. vanish in the limit as 
e— ~~ » whereas the Te are independent of € .) 


One might want to decouple each hierarchy as follows: 
(A), we (A) 
G + cam h. G,. (5-9) 
of “T= o 


where 


‘ ro 
L ae h independent of (5-9a) 


oh) 
Each set (labelled by 2) of (5-8) is now closed, and can be solved for és 
Then the sum (5-6) can be taken over all the solutions. The above procedure 
is not equivalent to that leading to equations (5-5), for if we sum equations 


(5-8) over A , we get the set of equations 
ui at, ~ G, 
e Oath toy 

; (5-10) 


Wee hoe 


Because of condition (5-9a), the set (5-10) is not closed, as was the set (5-5). 


e Og 


This method of decoupling is also unsatisfactory in most cases. What 


() 
usually happens is that on is also made up of elementary Green's functions 


C (AYA) 


«) 


G, Z.. mie (5-11) 
jf? 
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A CAB) ; 
and one might want to decompose om in terms of the ©, with different 
(3) 
coefficients for each @)(4) pair. So separate equations for os c are needed. 
Continuing on in this way, we find that our most general set of equations can 


be written 


A) scene” CAD ieee (A) 
2G, = “+ = G, ‘ 
fs (5-12a,) 
Ae) 4) CNAYED 


e G, = ee Gs, (5-12a,) 


(A)... (8) Wd ()..(S) « W),., De 
(A Wp a ae ” 


a 5-12 
ely EY s sa 
é 3 
ree each Green's function that appears is elementary. Here 

cn (MBM 

et » the result of some specified decomposition 
(VA) 1 (D) (€) oo 1 OD) 
is aed it (5-12b) 


is some linear combination of all the other elementary Green's functions. 
Equation (5-12a) represents UL equations, one for each value of & ; 
equation (5-12a,) represents ( JL Y¥U ) equations, where YVL is the number 
of (a -values; and so on. The number of equations in the set (5-12) is equal 
to the total number n of elementary Green's functions that appear. We call 
equation (5-12) our "closed set of decoupled equations". 

For convenience, we relabel all the elementary Green's functions in (5-12) 


by the single subscript eo . Equation (5-12) is then rewritten 


n 4 
eG, 2 Fy ok 2 ae G, API 2,u,N (5-13) 
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where Cow are coefficients which may be zero, and are independent of E. The 
Tax are the thermal averages associated with the Oy » in accordance with 
equations (2-50). 

In order to identify G, , we call S. the set of all values of 


such that 


G, = o Gr, (5-14) 


, TE, 8, 9 
in terms of the Gy and the G, . For, summing (5-13) over xé S, , we get 


aS Ge Ge + S00 es (5-15) 


eS. Aé Sz KE So 


It is then possible to identify Te a, Ti QO. 


where the summation convention (repeated indices are summed over all possible 
values) has been used, and will be used throughout the rest of the chapter. 
Comparing with (5-2a) for p=0, we get 


aos 2 


ALS: (5-16) 


G,  < Ge (5-17) 


Xd & Sp 
We also have, still using equations (5-13), 


(ji = J B Ai pFe 
© BO Gute BEC Te TBE Ga toa 


LES, ( 


Comparison with (5-2a) for p=1 gives 


(5-19) 


G. = ina Cy Or ; Ge, (5-20) 


4] 


In general, for (2g we can write 


Te a o0% er was 
Ae Pe (5-21) 
fe 
weak an ae “i Cra (5-22) 


as up to this point the decomposition (5-12b), which represents our (qt+l)st 
state decoupling, does not affect the original hierarchy of equations (5-2a). 
For {>g » the above identifications ee 21) and (5-22)) cannot be made. 
However, ie primed quantities — and a by 


LOW, ‘a fe fe (5-23) 
‘hs — fr Aw [PA , 
Ge = Sa Sa a Ges Ge 


(5-24) 


for all é » we can continue equations (5-15) and (5-18) into the infinite set 


/ ay, / 
a G Lge py , all = 


which is based on the closed set of decoupled equations (5-13). Note that 


ale: TT] «ks 9) (5-26a) 


and 


OG, (hed) 


(5-26b) 
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It is of interest to note that one can equivalently define Te ane G 
by equations (5-25), (5-26a), (5-26b), and the decoupled set of equations in 
which one is interested, e.g. (5-5) or (5-13). Equations (5-23) and (5-24), 
our previous definitons of Ty and oh » may then be deduced through the 
use of equations (5-13), (5-21), and (5-22). 

Let us consider, as an example of this alternate definition, the special 


case of decoupled equations (5-5). Looking at equation (5-5a,41)> we define 
4 


Cra by 
Gus, * ca Og 
1 420 
and write a Green's function equation for it, using equations (5-5) to do so: 
/ =| 
e © wee 2 dl, Ge +ed, G, 
Y £=0 4 | 
-| 4-! 
= S d, TV - al. an 
kre ~. 


a dT. a d, as oly G4 
peo 


* = d, Tp x (dp, +d, dy )Ge 


470 


where in (5-28) we define dias d_5».to be identically zero. We then define 


Wea by 
/ a 
leu = 3 aly lg (5-29) 


(5-27) 


(5-28) 


Og+a. = a (oh, +a Gy (5-30) 
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’ 
We next write an equation for Orgs, to get, in the same way as above, 


7 
Cc Cava . Fler ANT + (Green's functions) (5-31) 


£ 
and define Tosa by 


Pe = = (de. + , de) Te (5-32) 


and so on. 
It is this latter definition of the ly that is useful in practical 


applications of the Moment Theorem, which we are now prepared to state. 


3. Statement and Proof of Moment Theorem 


/ 
Theorem: The Q th moment Me associated with the line shape function 
derived from the closed set of decoupled equations (5-13) is given by 


/ 


Me Ay, [To (5-33) 


for all ¢ 
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Corollary I 


"4 4 
For (q+1)st-stage decoupling, the first q+l central moments Mos My a 
& 
Mare exact. 
q € 


The proof follows immediately from the definition (5-23) of /¢ and 
the identification (5-21) for P< 4 ; 


Corollary II 

If one wants a line shape function which gives the first qt+1 moments 
exactly, one cannot afford to decouple before the (q+1)8* stage. 

The proof is analogous to that for Corollary I, with the additional remark 


that a thermal average cannot in general be expressed exactly in terms of lower- 
order thermal averages. 


Proof of Theerem 

The proof will be given in part for a simple case first, and then in full 
for the general case. 

Case (a): Closed set of equations given by (5-5). 

Equations (5-5) form a special case of the set of equations (5-13). For 
proving this case it is more instructive to use the language of the special 
case, which means using the definitions of T; and Gr. given at the end 
of Section 2. 

Since equations (5-5) represent q+l linear equations, each , (o<t <4) 
can be written, by Cramer's Rule, as a quotient of determinants, the denominator 
being a polynomial one degree higher in G@ than the numerator. Barring 


accidental degeneracies, then, each Gy can be written as 


vi 


ai 
G, = ne en (5-34) 


Pl = Ep 
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which results from the decomposition into partial fractions of the quotient 
solution. af are certain coefficients independent of E, some of which 
may be zero. The at are rea] as is shown in Appendix III for the general 
case. The Ep are real, distinct, and independent of e . The reality 

of the Ep follows from the analyticity conditions on the Green's functions. 
Their independence of € follows from the observation that all the Gy have 
the same denominators in the quotient of determinants. 


Substitution of equation (5-34) into equation (5-5) yields 


r 
x EEA = T+ Aa gly 


r 
P 
(2-£,t€,)Ag - A 
Se SEE TY AE A. 


=| 
Equating coefficients of eé- Ep) results in three conditional equations: 


SA, = Te o < f$4 


(5-37) 


P f 
E, Ay = B ou ;_; = = {<¢ (5-38) 


_ P (5-39) 


Q=0 


Now, using equations (5-3), (5-34), (5-4a) and the asymptotic formula 


bin ho = PE) F ir) 


+! (5-40) 
e->ot MEve 


where P denotes the Cauchy principal value, we find 


* 
‘(e92-s) wry = 


(€p-S) / 
(€-S) woas Ay yg en er a 
mS lav d 
(2~-S) iv = fa ca. se P= eo ee) 
‘ % Sirs 9 — St *(6E-G) pue (ge-G) Wor 
; oJ 
ij 
(Lp-S) J a 
JX > 
Cz) re 
, d 
I3-to-1) GV < oy 
d e ° 2 - 
(13-01-08 IY LZ Corp cop 
3 =Hiae® “3- °c -3 40e3 
———— ee er li pl * sal a : 
ay ete a 13-7 -2 


fe + wi Wale ™) cpp “a o 

ome ee as doe 
40<3 = 

Brae aiig Tr Coe) cpl a 


O96 Corre ops = 


9” 
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So, for TS rf » equation (5-41) yields 


/ ) / 
Mp = lar / (5-44) 


If T= qtl, we have 


Sb PA.= oo E> AG from (5-38) 


_ =A Ve from (5-39) 


€20 
af > AT; from (5-37) 
€=0 
Pe a rol from (5-29) (5-45) 


q 


and so 


Z 7 / 
—— (5-46) 
M = | ce 7 / io) 
from (5-41), (5-44), and (5-45). The theorem is thus verified for 4 ae ee 


q+1 (equation (5-33)). The proof is not continued for higher values of 7 for 


this special case, as the notation becomes quite unwieldy. 


Case (b): General closed set of equations (5-13). 


_— 
For this case we make use of the definitions (5-23) and (5-24) of 7g and 
/ 


G, « 
From the theory of linear equations, the solution of equation (5-13) is 
given as a quotient of determinants, the denominator being a polynomial one 
degree higher in @ than the numerator. Barring accidental degeneracies then, 


each G, can be written as 


- 
nD OA 
eet are ae (5-47) 
pa 2 
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which results from the decomposition into partial fractions of the quotient 
solution. n, the number of elementary Green's functions, is the degree of 
the polynomial comprising the denominator; rv are certain coefficients 
independent of E , some of which may be zero, and the Ep are real, distinct, 
and independent of « . The reality of the Ep follows from the condition 
that the Green's functions be analytic on the upper and lower half-planes. 
Their independence of follows from the observation that all the Om 
have the same denominators in the above-mentioned quotient of determinants. 


Substitution of (5-47) into (5-13) yields 


P 
(4- €,te, ) A we a he + —o As (5-48) 


=| 
where the summation convention is still in use. Equating coefficients of (e-€,) 


results in two conditional equations: 
PS 
) ie * Bl} 0% (5-49) 
fr 
Zw ve ch at (5-50) 
‘aa 4 Ai 
A at | 
It is important for the proof that the ol be shown to be real. This 
is done in Appendix IV. 


a 
Multiplying equation (5-50) by 1, Ep 1 Head eiens successively, and using 
(5-49), one finds 
a 


' A 
2 & A Ae COE Ags CLT, 
= 6TAls Be Sih = CAG 
Po, P 


- 
Seat Ce i Oe ae, 5-51 


From equations (5-3), (5-4a), (5-4b) for p=0, (5-47), (5-40), (5-14), 
and the reality of the f » the moments of * (i) are given by 
BA 


Me = bis, (do Guy oe, OE Pps sae it 
Lime (a) Tm Gyleers-te 
Lim 3 Ean re ©; 
e—>» oF [os Sh, aL a A lesw-le 


Liw 


E-9 o* (alo tn 2 Galesu-it 
“im Sd (os~w2)* = = TAs /e-14-£,) Bruais, 
Lim - fav a = (Ad /e-W,-8)) leew-ie 


fd (0. EZ >. Al S (y=, = £,) 
fl SS AL Pla G) 


AE Se 


» , ‘a ae by (5-51) 
AE Ss a pee aa 


by definition (5-29) 


This completes our proof. 


4. Summary on Use of Theorem. 


It may be worthwhile here to summarize this chapter indicating how the 
theorem would be used. Examples of its use will be given in the following 
chapter. 

The line shape function can be expressed in terms of a Green's function 
G, and the equation for G, leads to a hierarchy of equations involving 
higher-order Green's functions Gs - This hierarchy is decoupled in some 
manner leading to a closed set of equations which can be solved for the 
approximate G, - One could then find the moments of the line shape by 
suitable integrations with the approximate line shape function. The theorem 
shows that this integration is not necessary. The approximate moments up to 
the stage q+l at which the decoupling was done are just the usual thermal 
averages appearing in the equations and agree with the exact moments. 

To find the higher moments of the approximate line shape one continues the 
closed hierarchy of equations. In this continuation no new Green's functions 
appear; each time a Gy,, Green's function appears in an equation it is 
decomposed with the rule first used. The thermal averages appearing in these 
higher equations are the appropriate higher moments of the approximate line 
shape function. 

This procedure for finding moments has the advantage that it does not 
require one to solve equation (5-13) explicitly; one need specify only the 


decomposition (5-12b) that is being used. 
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CHAPTER VI 


Resonance Line. Shape Problem 


i Fe Decoupling Methods 


At the beginning of Chapter IV, the properties demanded of a suitable 
decoupling technique are outlined. These are the tractability, degree of 
complexity, and faithfulness of the resulting closed set of linear equations. 
It is, in any specific case, a simple matter to determine whether a set of 
equations is complex or tractable. It is the question of faithfulness that 
will be dealt with in this chapter. 

We postulate a faithfulness criterion: one decoupling procedure is said 
to be more faithful than another if the first few moments of its consequent 
line shape are closer to the true moments in some valid approximation. 

The phrase "in some valid approximation" is inserted, as expressions 
for the moments are invariably given in terms of thermal averages; and even 
the thermal average of cis has never been given exactly. 

The criterion is practical, for the moment theorem of the last chapter 
immediately enables one to write down expressions for the moments which are 
to be compared with the true moments. In line with this practical approach, 


only the first few moments are usually asked for, as higher moments become 
increasingly difficult to calculate in any approximation. 


Our criterion will now be used to investigate the relative faithfulness 
of two specific decoupling methods, when these are applied to the problem 
of finding the paramagnetic resonance line shape function associated with 
the physical system described by the no-crystalline field Hamiltonian (2-9), 
where “WP and Wy are given by equations (2-11) and (2-12a) respectively. 
One method we shall call the excess-over-random (EOR) method; the other, due 
to Tomita and Tanaka (1963), which is a partial excess-over-random.method, we shall 


refer to as the T-T method. 
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We specialize our physical system still further by requiring that it 


consist of spins having S = 4%. For this case, from 


DORE ot Star (6-1) 


and the well-known commutation relations for spin operators, the following 
relations can be proved: 


GS ac SS et oS E 
So ores =n ee = ~ 2s” 
65-4 = - SiS - " (6-2) 


SSE a SOS 


The significance of the above relations is that any product of two spin 
operators at one lattice site can be reduced to a linear combination: of 


single spin. operators. 


The EOR method, which takes full advantage of relations (6-2) 


» proceeds 
as follows: 


One rewrites the Green's function equations (2-50) for the 


"ordinary" Green's functions [Al = [a/b] + € A/B 
comulant Green's functions (A/ = 


Cl = Gl 


in terms of 


(A/B) which are defined by 


it 


(ial = <p (al + <arcil + Ci2| 
fiaz|)= <7 (22}+ @isl + <39 (i2| Py 
+ 237 (1] + isr(2al HKD + (23 


etc., where the numerals "1", "2",... represent various operators. After 


oo 


the above decompositions, one stops at a certain order and discards the 
highest-order cumulant Green's function. For example, in second-stage 
decoupling, one stops after equation (6-3c) and neglects the cumulant Green's 
function (1231, supposing it to be small. One of course verifies that the 
operators "1" and "2" in, for instance, C123| » are associated with different 
lattice sites, Otherwise, the product of operators "1" and "2" is to be 
treated as a single operator "12" in the above decompositions. As an example, 
the Green's function (5°5°5" 1 (ia J) for the case S = %, must 

be changed to va £5*| before any decompositions are made. 


The T-T method proceeds in the same manner, with the one essential 
difference that products of operators associated with the same lattice 


site are not treated as a single product operator in the decompositions (6-3). 
The first two Green's function equations, written for the elementary 


Green's functions, are 


€ KIER He CSD EN 6 PIED, 
ECR SCR OM 
2K LIER = BEL ED* MSGS ISR 
EO €GSG SR - FEA ES SSE 
++ > Ae GHG IER  (6-4b) 
Note that we in fact want to Find 


G2 BLE ISR 


(6-5) 


These equations hold for any effective spin S. 
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Second-stage decoupling by the EOR method, for the case S=1/2, consists 


in making the following replacements: 
kG 5G SEM, Garay Pp CI GER | 
th G1ER* PED -B)EGGN 
ES GS SoM ale], $577 5" 1S Xe (6-6b) 
+ LOGIE G EN 


~ 
where J = CS: 7 is independent of i by translational invariance. 


6-6a) 


When any two of p,i, and j are equal, use of equations (6-2) ensures that 
no new Green's functions are present, and hence that the resulting set of 
decoupled equations is closed. 
Second-stage decoupling by the T-T method consists in making the 
replacements 
aorertic- Rertic- 
CSP S*5*1GR aaa a ¢€f, Sj AF » th €S S or 
6-7 
+ (<g> - 2-Yqsism | 
CSP <7 ES LeMe 
- Ta C= (6-7b 
+ GD EG NR : 


for all p,j, and i. 

It may be seen that in the EOR method certain correlations, namely auto- 
correlations, are retained in their exact form, while in the T-T method they 
are approximated. Thus, if we take p = i#j in equation (6-7a), the 


replacement 


LGC > an ESIC R+ EPA) ET IGM (6-10) 
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is effectively performed in the T-T method; in the EOR method, the left-hand 
side of (6-7c) is retained without approximation. 


Substituting equations (6-6) and (6-7) into equation (6-4), the EOR- 


decoupled equations are 


e[tl= beg + Ay [B71 = cy [FT| 


(6-8a) 


and 


e[Fil= oe (ad SSP <I <Q) + 7A) Li 
wt oa RUF T-4Cy [tl + 4A. re 
4 SEAM te ZA, Dy CPT 
a (\- 0} ip y Walt Paksl +(S, tty ayn}tif 
26 ORDA Tels TR )-3- tf 
“SEA, (1-8) - Be) {<GGOUEL EKG OL 
. <5 at 7 SG Str i/ (6-8b) 


where we have used the notation 


[tl = CST 
pene <G SS. 


(6-9) 


e[tl-+ Rey + ZAG EZ <j [Ft (6-10a) 
© [Ete de de <tt> ~ Fe CG") 
tEAM Sal? 2 Hi [tA] ?, t+ ((s26%>-a.-)/tf 
Ee, Lelethpls (GEO It} 
. LEA Ls SG [t+ Gs, vats 
Gt FGF 


(6-10b) 


Other types of second-stage decoupling may certainly be tried. One type 
that seems worthy of consideration is due to Callen (1963). It makes use of 
the fact that =" may be expressed in two ways in terms of as and 
ul > “for Sie kb 18; 


=. 2/fotrrs. pee 
a (S f i; S, ) (6-11a) 
— Per Ge eee (6-11b) 
‘ : rc : ‘ ‘ 
f. "hs replaced by 


Era. fe (SS- 676 f\+ (l-u) (t- £7) (6-11c) 
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after which decoupling proceeds by the usual EOR method. Here o{ is an 
arbitrary parameter which must be fixed by some suitable physical criterion. 
The Callen method will not be considered in this thesis, as the calculations 
become very long and involved. 

Other possible methods have been looked at and tested on the Ising 
model Hamiltonian, but they were found to be of far less value than the 
methods being considered here. 

From Corollary I to the Moment Theorem, the first moment, derived from 
both the EOR- and the T-T-decoupled equations (6-8) and (6-10), is exact, 
since q = 1. Let us look,then, at the second moment. By use of the procedure 


Suggested by the Moment Theorem, one arrives at the following expressions 


for M,(T-T) and M,(EOR) : 


M(-T) = - oF oe sa Ef BABy (G2 +S HY 


+ eG 5) = why (Kot +Ke5*)) 
—p Ay (KEG P+ SE SI- GBP-K GLY) 
te Cu (KESIP+ KGS) +4 Bi} 


(6-12) 


ay eee en se Flu) {<G Ss A (6-13) 
He Rly LG RON so rapt 


The correct expression for Mo» denoted here by M,(CORRECT) , is derived from 


equation (4-15) with 


bon Se a (6-14a) 
B= a (6-14b) 


It is given by 


My(CORRECT) = — oh Seba (a LS P57 
~<4*s- ‘5°? CSESPS) = Cy ASL GE? 
— ESS >- eG" *) - Ag keg > 
+e Si Ps 7 = SF fe S: WA rede) 
- <Se*S"7)} 


To compare these formulae for the second central moment, we use the 
McMillan and Opechowski approximation technique which was described in 
Chapter IV in the section after equation (4-35). The following relations, 


valid to (MO), and for S = %, are needed for our proposed comparison. 
fs pe thE (par id) | (a= i/kar) 
>> “= Fact Yo - 14)” 
= ae masini - /4) 

7m 4a Gu d)[@E-Z)h “hz 


(6-16) 
65") Axe yy MW 1a) 


wo . 9) 
where ws is 57 feh™ 6°) /T ef™ 


is the first-approximation part of , and is independent of i, 
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Gq =a. jie 
€ 


is independent of i by translational invariance, 


and Ze + Coe + © 


Substituting the above expressions into (6-12), (6-13), and (6-15) the 


expressions for M, become, in (MO), , 
2 1 


M..(e£oR) =M, (coareer) 
sae (ay! Baty + tel YS yf, 


M,, (T-T) 


J (6-17) 


in gts 


(Tat) = RMS of (619) = soe ye W) Fe Buty Lainey WCig 
+ Cy jt Ge)+ + = (Rie + Ai) .* 
: oes 


and 


Mo(EOR) = Mp(CORRECT) = R.H.S. of (6-18) + 2/3 (Me = 4) x 


am TS 
Si, 
‘ 2 (Bd fons Halie es, (6-19) 
t 
And so the second-approximation part of M,(T-T ) is found to be in 


error. Moreover, this error is not 0(1/N) of the correct second-approximation 


part. This may be seen by comparing, for example, the terms 


and 
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The first term is not 0(1/N) of the second even though the latter involves 
an extra summation. This is so because Big and Cy¢ fall off as rapidly 
as V/in-2}? with increasing distance between lattice sites. 
O'Reilly and Tsang (1962) have done summations like the above for the 
special case of Calcium Fluoride crystals, bearing out this statement. 


Note, however, the error.is nil when T = 0, as in this limit 


and also when T = © , as then 


peer 


We now look at third central moments. With considerable work it can 
be shown that, to (MO), 


(EOR) = oti 


= = fam," @) - WV. Uke? 4) [ ABB +2 6 


(6-20) 
-33 8 -A BB +77 a A Bi} 


ie 


where the notation 


ABCD... - SH Ag Bele Okun 
‘4 


has been used. 
M3(T-T) is more difficult to compute, even to (MO), . but it can be shown to 


be in error generatly by showing it to be so for the Ising model. For this 


case 

7 2 pn 
M, (CORRECT) = eee 4M € ee fr] (6-21) 
and 


3 v 
M,(T-T) = M,(CORRECT) - bé N Me (A. ' /4) (6-22) 
arT, 
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Again it is seen that the error is not of 0(1/N), and that it vanishes only 
in the limits T=0 and T = %, 
The following table is drawn up in order to highlight the differences in 


moments calculated by the EOR and T-T methods. The ticks represent areas 


of agreement with the correct values, the crosses, disagreement. 


By our criterion, it follows that EOR second-stage decoupling should yield 


a more faithful line shape function than the method of Tomita and Tanaka. 


2. Solution of Equations. 


When one comes to consider the solution-of the EOR and T-T-decoupled 
equations (6-8) and (6-10), the tractability requirement assumes extreme 
importance. The T-T equations turn out to be completely admissible of 
an analytic solution in Fourier-transformed space, while the EOR equations, 
which are almost identical, have some added terms which increase the 
number of simultaneous equations which must be solved for G, to 
something of the order of N. 

Moreover, the T-T method yields a set of tractable equations for spin 
values higher than 43, and can easily be extended to third-stage decoupling. 
The EOR method, on the other hand, becomes prohibitively complicated when 


higher spin values and higher-stage decoupling are used. 
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For the reason that T-T second-stage decoupling does not conform 
sufficiently to the complexity requirement at high temperatures(i.e., 
temperatures in the paramagnetic range), Tomita and Tanaka extended their 
calculations to include third-stage decoupling. Now, T-T third-stage 
decoupling leads to equations which,: according to our criterion, are 
even more faithful than EOR second-stage decoupling. This follows 
from Corollary I to the moment theorem, which tells us that the Mo» as 
well as the Mi» resulting from a third-stage decoupling, is exact. For the 
above reasons, no serious attempt will be made to solve the second-stage 
EOR-decoupled equations other than formally. A simple case will be 
considered, however, to show that these equations satisfy the complexity 
requirement, even at high temperatures. 

The first two Fourier-transformed second-stage I-T-decoupled equations, 


derived from the spatial equations (6-10) are 


fe - Bg (ey}(a| = + io D(g, h-4) (oa4 4 (6-23a) 
and 

fe - A, G)} (hcg Me Dag 4 +5 R(k,g) (t | (6-23b) 

where Do (4)= Wot <e? D@o)= Wht pd (¢,0) (6-24a) 

Rk.g= D(k,g-k) blag gel + CG REG F nae 

Digi? Gi 6G) (6-24c) 

C44 = a (AG)- A @)) (6-24d) 


iba. Cie 
A (k) = Wa 2 & Ay (6-24e) 


and the notation of Tomita and Tanaka 


Cel = (tlaj= fel! 
= a= eth {f oT al: Zee ey 1S yy (6-248) 


(keg. (3 ‘k) 
ea 


($3 Iz) E the)- seria! é} (6-24h) 
<5 See - “<8? C51 SER 


has been employed. 


(6-24g) 


Li ‘Das q Veg 4- Hh and tag are Fourier transforms of 
certain linear eile of thermal averages. 
« /H= ft - 
Note that G, (t= (flee (6-25) 
Equations (6-23) are easily solved for (t| . The solution is 
+ a+ Si) 
(t| a 4 E- be (9) (6-26) 
a E- O,(k) - Salk, 4) E -b.(9)] 
where ¥ 
= , 
a, (le, 4) = y Dts, k-9) Aung, 9 (6-27a) 
is ge, a 7 
Q (kg) = D (4, ke-g) RC, 4) | (6-27b) 


Equation (6-26) does not satisfy the complexity requirement for high 


temperatures, as 
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bs. (4) amy ty), ‘i barge) 


(6-28) 


and an (vo) reduces to a sum of two delta-functions. 
The first two Fourier-transformed second-stage EQR-decoupled equations, 


derived from the spatial equations (6-8) are 


fe-a,ui(t]* 4, + (ED (5, Nee i (6-29a) 
(e- OF big g |= 4443 + UR ei gNCEI 


+ [ln PRN SGM, EI] co-2m 


where 
Plkgi= Yn) 5 1C Ceekg, hegeh’) Ze (6-200 
+ Clherh’-g R) thy r= D hgh’) Ble ee} 
Qlkgk)= ED Ck’ k-k’)= 2% DEK k'-g’ 
<A>Dllrkng, R-k’) (6-30 


To solve equations (6-29) formally for eq » we treat the terms in 


) 


the square bracket of equation (6-29b) as sufficiently small, compared to 

+ ; , or 
(WN) R (kg) (ft | so that an iteration procedure can be used; this 
is equivalent to the recognition that the EOR and T-T methods do not diverge 


widely from one another. Letting 


Rlk,g)= Rlleg) + Plk,9) 


) (6-31) 


the second EOR equation, (6-29b), may be written 
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te-S@} (cy 4] = Qing g + WR OAG) AI 
+ (1/n) SS high) A pl (6-32) 


i 
or, with the use of the iteration procedure on (eo 9 


\e- A. (9) (24 |= S? (k,4)+ (An) fe (ha) ( (F | (6-33) 


where 
= ‘ 4 \/N Qlle 2D (le, 4, h 

Se lkigy= Shg.g + (a) Se a ) dew t 

+ (\/n*) > Qlk h) (kh, h") Shek pi h! 
Rok" = (e-0,h)(E- 0,(k) 
(6-34) 
ao tag 

We note that Je (k 4) satisfies the "integral equation" 

S, (be, h’) 
Sz (e.g) = Ang, 4 , + (ANE Bas — at) (6-35) 


Similarly, oe (4) satisfies the "integral equation" 
Clea) sf (k, k’) 
‘ came ESS ) ' (6-36) 
Jelhghe Ries) + OMI Ty 
Equations (6-29), solved simultaneously, yield the formal expression 
+ we D(g,h- 2) & & Oe. 5) 


E- 4, (4) 
— en es ee . 
( = E-A, Oe DG,k- 4) de (k 4) (6-37) 
E- D. (3) 
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One would like to solve equations (6-35) and (6-36) exactly for 
ie Cie 4) and of (kh 4) respectively. Unfortunately, this 
‘ 
cannot readily be done, as the coefficients in these equations do not 


possess convenient symmetry properties. Let us, however, consider the simple 


case when 
Qlk ah) = Qike.a4) dhe ) (6-38) 
Equations (6-35) and (6-36) then have the solutions 
ne 
le io Fie, 2 ass oe 


R (kh 4) aft 
J (k.4) nt soa = ae 


E |= 


Substituting (6-39) and (6-40) into equation (6-37), we get the following 
expression for (rt 


Oe +y Dg keg) ding g 
N E- 4,4) = = Q (h58) 
PO A a Se ) oe ee 
k te lain al D (4, h-g) R (k, §) 
N* Fo Es Gol) Olk,5,8) 
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Equation (6-41) satisfies the complexity requirement even in the high-temperature 
limit for which (6-28) holds. We note that this requirement is also satisfied 
when equation (6-41) is specialized to the Ising model case, when (see equation 


(6-47) ) 
4A, (9) = 8, independent of q (6-42) 


The latter statement cannot be made about equation (6-26), the 
equivalent T-T solution. It will be demonstrated, in the next section, that 
also 3rd-stage T-T decoupling is not "complex" enough for the Ising model 


case. 


3. Complexity Requi-rement 

We have seen above that at second-stage decoupling the EOR method 
satisfies the complexity requirement that the line shape consist of many 
delta-functions. In order to have this requirement satisfied with Tomita- 
Tanaka decoupling, these authors had to go to a third-stage decoupling. In 
this section, we shall discuss this point. 

Using third-stage decoupling, Tomita and Tanaka arrived at a closed 
set of equations which proved to be quite tractable, having as a solution 


kl = ur /pr 


(6-43) 


with 
o *(h ¢ 4’) 
+ av (kg) — 2 = = een Te 
Jy = E- Aglk\ - > qé = Ba) , (6-44a) 


E- 4, (3\- b*(h, 4,39 
% ‘ 7 &- b, (¢”) 
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* E-o E-0 (9) 
E-o.4)- 5 bothag? 


pe E- A, (4) 


we at iis co" (kg) 
ane es 


(6-44b) 


where Co (k,9.9’) and co (k, 4,4’) involve Fourier transforms of Aig» 


Cig and of certain linear combinations of thermal averages, and 


b* (k,4,4') = D(4',9-3’) 3D (4,459) dg 544 “CQ a 
+ 2G’ k-g- 4) ] D(q, bg) 44 


~D(qrgth k-g) Vegeg? eg th . (6445) 


Equation (6-43) can, by use of the method of partial fractions, be 
cast into a form which leads to a sum of delta-functions for the line shape. 
These must be smeared in some way so that a continuous line shape function 
results. One can perform this smearing only if the delta-functions are 
sufficiently large in number, and their arguments sufficiently dense 
along the W-axis. It is this condition that has been referred to as 
the complexity requirement. # 

At first sight, Tomita and Tanaka's expression (6-43) for (Ti seems 
to satisfy the complexity requirement. It especially seems to do so for the 
general paramagnetic case. However, it demonstrably does not for the case 
of the general Ising model, at any temperature and for any spin. This will 
now be shown. 

# Grant and Strandberg (1964) satisfied this complexity requirement by 


treating the distribution of spins well-removed from the one of 
interest as continuous. 
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Recall that for the general Ising model (2-27) 


A. = 
| ¢ (6-46a) 


A (4\5 O all (6-46b) 
As a consequence, the quantities 4, (g) and “Ch 4.40 involved in 


equation (6-43) become 


A, (¢)= Ww, + BoM A, independent of q (6-47) 
and 
bk,4.4))= ‘a (4-9") Ving arg 3 b (ke, 4-g’) (6-48) 


Then DT and oa (equations (6-44a) and (6-44b) ) become 


3 a*(k,4)- (E-a.)' 5c (k,4') 
4 cote 


+ at O,(k,g) + (E-a,)' SC, Ck, 4,4’) 
be 5 ¥ <a j 
E-A,- (E- a Sb (k, 3-3 ) 


(6-49b) 


22 b (le cop a a b (k= b* (RY (6-50) 


3 
by translational invariance. So, if we define Gla) Cok) val (k) and 
Qn (k) as follows 


70 
ad) = z— HER §) 
OPO Se ee SO AY 
—$ 4’ 
Qa (kh) = = a*(h, 9) 


(6-51) 
ots SF 60% 5,3) 
we may write 
gt 4. Gethy + (E- ayle.*(h) 
.. ; E-o, - (E-a,)"! b*Cke 
( ‘Tah ae Ea 
E-a, - (€-a.y! bYlh) (6-52) 


This yields, via (6-25 ) and (5-40 ), a line shape function which is expressed 
as a sum of three delta-functions. It is, of course, wrong to extrapolate 
from these to a continuous line shape function. 

Now, the general Ising model admits interaction of any spin with all 
its neighbours, no matter how far away. The most elementary considerations 
indicate that unless one is dealing with the one-dimensional Ising model with 


S = 4, nearest-neighbour interactions, many more than three delta functions 


7| 


are required to describe the line shape. One must conclude that Tomita 
and Tanaka's third-stage decoupling cannot be used on the Ising model with 


any success. 


4, Application to Case of Completely Non-equidistant Single-Spin 
Energy Levels. 


This chapter so far has considered the case in which the single-spin 
energy levels are equally spaced. The most fruitful approach to the case 


of completely non-equidistant energy levels would appear to be through 
the use of the operators ai, introduced in Chapter II. 

This, however, proves sufficiently complicated, even with T-T decoupling, 
that it is unlikely to yield a line shape, and one has to be content with 
calculations of moments as in Chapter IV. To see this, let us discuss 
first-stage decoupling and indicate, what would happen in second-stage 
decoupling. 

The statement that the unperturbed energy levels are gompletely non- 


equidistant implies that whenever 


one has 


A” Sag Gh (6-53) 


This is a restriction on the elements of the set G defined in Chapter II. 


The Hamiltonian my can then be written as (see equation (2-23)) 
_ a “ + 2pzt 
RecN, * SE} (a5 4%, a FG Gf (6-54) 
U eal A“ 7 d 


(0) 
where ‘MH, is defined by equation (2-16). 


t 
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We will take the unperturbed frequency, W, ,of the primary line under 


. consideration to be given by 


bly a gt, (6-55) 


We than look for 
Ge hones d_ oe CaN Agha ye (6-56) 


in order to find the line shape function 


(6-57) 


4 (3) = E> ot 
cin > Jd sae Oo lumte 


g—>o 


The first Green's function equation is 


e (ayia Lan, wal 


ey Amt of Lada Sanda 
d 


+ ae ee Bir. allel 
= = Ar 4; Adin. Auin)at Ay Aen Aun Nel An) | 
d 


- * 
x =o; 2 Cm [An lay, 4jn]| uae) 
d 


where we have used the relation 
z > = 
$2 EEA, ai) 


= (6-59) 
= a A: } 
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which follows from 
“8 - = 
[ s S-] = 2S: 
2 
combined with (2-32) and (2-22). 
T-T decoupling in the first stage consists, here, in performing the 


decomposition (6-3c) and neglecting higher-order cumulant Green's functions. 


This is equivalent to making replacements like the following: 
caw ers. |— Cait Anre) (4;4| 
- + 
- ats. 8 A) (4.4 | ) 
[Anis oe Aj; = CAmdei amr (Aji | 
© CAinres Afi? (Sand | 
With this decoupling, equation (6-58) becomes 
€ (at, | = — 4 [Aik Agia) * ZA 8 


+ = = a <4) aes. (2 a 


- — a C4j\ea Amar (ae | 
” 2 Cm coy (wn | 
d 


(6-60) 
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Taking Fourier transforms as in Section 2 and using (6-56) and (6-57), one 


gets a single delta-function for the line Shape: 
hays { = WW, ~ al C (0) 
2 y WW © fe A (k) tA 


- Wc . (Chi Waa) ~ Gea ies lf 
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where Via ‘ : Ayy independent of i 
pp ‘ a independent of i 
ipa 4 _ ae ~ ths - + 
Cin) Kay ” 2.28 c C dink jiu) : 
— 


In deriving equation (6-61), one had simply to solve one equation for 
‘ / 
feo ¥ bhp + 
Cir | = 2 & (4,, 
d 


To see the difficulties that arise when second-stage decoupling is 


attempted, consider the simplest case, that of effective spin S = 1 (with 


A ety say, in (6-55)). One must keep all Green's functions of the form 


- + 
' Zam, oi rad 


) 


—4n Ape, 


and one 4 in each Green's function. There 


si + on 
Avs yi Bag 
ane, pe: ’ 

where there are two 4. ’S 
are nine of these, seven of which are independent. A separate Green's function 
equation must be written for each of these seven, a decomposition made, on 
the right hand side of each equation, of higher-order Green's functions into 
lower-order ones, according to equations (6-3), and Fourier transforms taken. 


One then has, optimistically, seven coupled simultaneous linear equations to 


75 


solve. The latter number increases quickly when higher spins are used. This 
process is quite unwieldy, and so the analysis has not been pursued any 


further. 


5. Summary of Chapter VI. 


We have, in this chapter, applied a moment theorem to compare the 
faithfulness of two decoupling techniques. The results of the comparison 
are summarized in the table on page 61 . The EOR decoupling proves more 
faithful than T-T decoupling, but unfortunately, less tractable. It 
satisfies the complexity requirement at second-stage decoupling, whereas 
T-T decoupling satisfies this requirement only with a third-stage decoupling. 
For the Ising model, this complexity is not satisfied even with 3rd-stage 
T-T decoupling. 

The argument in favour of Tomita and Tanaka's approach to the general 
paramagnetic resonance line shape problem is that it does, at least in the 
case of equally-spaced unperturbed energy levels, provide an analytic line 
shape function. This may be more important than its providing correct 


moments. 


APPENDIX I 


To show that in 


CCS ‘ e i —_ Ata e ui rate (A-1) 
NEG, ) 


Se Ast may be replaced by SD Ary4, , 
d = ‘a NGG, “i 


Now, 


(ret). 2 e der Ab, reeks Balt 


(A-2) 


where we have expanded the trace in terms of eigenstates lL? of 7 
with eigenvalues Eu . 
Let 


jb oD 2itit 
Mi = ult Ajxne "2 pa Py, (A-3) 


Then 


OE | 
E mre eet Aine Ain ivr W| y, (A-4) 


Using (2- ea (2-21), and the hen eich ty of 1 » we have 
Ey Myy = <4|W | “Am, ee, i 2 EA? 


+ dd fet as ‘oe POT of WE 9465] [7 
th 
= 2. wit +QWer Manne i Seat Hye me 
+. dale? ieee yes Ai eit Ik? 


SEM + ToT — (tT ) Mi 


7 
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Therefore 
nl Se ag. aye 
a (Th . )~ (Trir rey Myr =O, (A-6) 
But, for AEG,» ee ae = WwW, (see equation (2-31)). 
So 
ak Xje 
tN gs ToT, My 
(\€G,) 7 as (A-7) 
Hu / 
18.5 My (Aé G,) is zero unless A é G,, 
Similarly, 
Nv ne 
XA Re, tte ye eA (A-8) 
(\€G,) a 
where 
AL ~~ Mot + iwor¢ 7 
Nv, = Ce c As c. Satz: EY, (A-9) 


From (A-8), (A-7), (A-9), (A-3), and (A-2), 


_ pps mn We + we % (A-10) 
rs <b Ae et ge abe 
A’ EG, ) AEG, 
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APPENDIX II 


Calculation of Thermal Averages for the One-Dimensional 
Ising Model using a Difference Equation Technique. 


As an illustration of the technique, the derivation of the thermal 


average 2 <'7_ will be given. 


2 
Sl) So = CS; 7? independent of i 
= TA (e-7* <n) where we have taken i = N 


= D> _ Mn). nan My) ss (Maer) at) 


nina <n Sty 


where n> Yi; are eigenstates and eigenvalues of or. 110 
) 


eager, = n/"> (A-12) 


and YI: can assume the values + Une j a = \ [eT ; and the 


following notation is used: 


n)= cath lay | is Ae ; (A-13) 


Pee ate Ae 
(1: Yh.) = Onli’ Se hi... > 3 CAC (A-14) 
We will also use 


+5,* 
at? a Nort ny (A-15a) 


iio 7 = <n, | — shy (A-15b) 
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f= ot bt (A16a), 
yer 4h (A-16b) 
f* at is (A-16c) 
Pe cthb Oy . (A-16d) 


Performing the summation over Vy , 


ed = a - 1). (Mn OHM su (Maer Myer) 


[a bn rcinty-£ b* ny ned , 


(A-17) 


Looking at the first term on the R.H.S. of (A-17), we write it in the 


following way: 


(2<1), = L* 2 <r rOnnn) MY) Ms) w, 
I P= Wei, You 
“ (Yivr) Mn Mn yor) Mtr 7. (A-18) 
Noting that 
Se XE? (Mi) 1 Minn) = Sat P+ rons? (A-19a) 
ni 


and 


= <UPOny(n, Mian) ot Oe +0? (A-19b) 


it follows that equation (A-18), after it is summed over ni, Mein Me (p ¢W-2) : 
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may be written in the form 


(2 <'9) - be po wae (Bia <uh? ” Bi <p>) 


Novi ee N yal 


’ (Meas) (“ine pvr) -— Eee, (A-20) 


where Gp and b, satisfy the recursion relations 


Q = X Gh, +/3 be 


(A-21a) 


a = {ant yb, (A-21b) 


with the initial conditions b, = 0, a, = ir 
Equations (A-21) can easily be cast into the form of a pair of second-order 


difference equations: 


ban (Mba = (Yn ad) bee = O with beo, bey , 


(A-22a) 
Apa ~ (ATA, = (Ha-AP) Op 29 vith Gel ara 
(A-22b) 
Equations (A-22) possess the solutions 
. p P 
a, = FL AY + G0.) mck 
Law Reset, (ay 
uate Fou big (A-24b) 


where Fo Fy Gas G. 


initial conditions, and A+ are given by equation (3-13). 


are independent of p and are determined from the 
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In the calculation of the second term of equation (A-17), one has 


occasion to look at the equations 


Cyt A Cy +/ a, (A-25a) 
don = YEpt Pd, | (A-25b) 


with initial conditions C, = 0, d, = ], 


They possess solutions 


ce FO eG. O.! (A-26a) 


dd. = Fy BS! = GA (.)F (A-26b) 


where Fo Gas F Gy are again determined from the initial conditions. 


d° 

Whatever thermal averages are done, one can always arrange the terms 
in the appropriate summations in such a way that either (A-21) or (A-25) 
need be used for each term. A list of the values of the F's and G's 


follows, where the notation 


Ye ) (4-0) + 430 (A-27) 


is used, 
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(a-Sar) MeV 2 eA 


ay Harter ~2Y¥ [(de d= PK 


a3 H+ o+V- ‘ = 3 / (44 PK 


(hear) [Ga i (aD /Gr AH 


From the above table, the following relations can easily be constructed and 


are found to be useful: 


dps 4p + a by (A-28a) 


(- b= $e (A-28b) 


Continuing from equation (A-20), performing the final summations, and 


using equations (A-21), we have 


(aem,= F (Gua <n? + By <iniay) (va) Sv Y 


a ~ One, — 


Similarly, the second term of (A-17), which we call (2 <7), ., is given by 


(2<9), = a lv (A-30) 


83 


and so 


Zk (= <'7), a (2<I>), 
Ve bain ~ Ane) 


- =: = by. fy 0m (A- 28a) , 


(A-31) 


One similarly finds that the partition function Z is given by 


Ze Gazi t Avi 
= (PO) + Gt Gu) (A) 
- eS Gan — 


as given by Huang (1963). 
All other thermal averages can be found in the same way, when this is 
2 2 
tractable. For example, the thermal average <S. dices becomes 
cumbersometo find by this method, as eleven summations must be done 


explicitly before the difference equation solution can be applied. 
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APPENDIX III 


Proof of (4-15) directly from (5-2a) and (5-3) 


Let i) n 
F, (+) = a ak G, () : (A-33) 
A 
Te * 2 Gey ie (A-34) 
Then = m 
(= Ge). (838) 


, u 
Multiplying both sides of equation (5-2a) by (it) || and summing over n 


from 0 to *_ , one derives the equation 


e F(t)= TH) -L mh Fe (t) (A-36) 


which has as solution we 
, iedt-f 
Fes i fa err | (A-37) 


A boundary condition is needed to fix to, the constant of integration. We 
make use of the fact that when “‘ = QO , the hierarchy of equations 


(5-2a) reduces to the single equation 


eG, 


r (A-38) 


and moreover, 


Tipe" T, 


‘ (A-39) 
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Substituting (A-39) into (A-37) and using (A-35), 

é: Gey s Ty (i-e 9 (A-40) 
when ”"” = 0 . To make (A-40) correspond to (A-38), one needs 
Ce+w (In & <0) (A-41) 


The assumption is then made that (A-41) holds also when yu 2-5 


Now 


, (ih 3 le “En Gitte 1 2 
RA E> ot 
im Low t ( 
4 ot Fonte 2) 


i's aft fal Rifles 
af Ate eT 


= Lim + fag ripjerel- c |f| 


2 ot 


_ a, £ (ap rer eeltl 


f» ot "ob he 1°] h! 


(A-42) 
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where we have used the relation 
ee 
e pent & f) (A-43) 


which follows from (A-34) and the reality of the fgg - Equation (A-42) 
is equivalent to the stage reached in Chapter IV, equation (4-7), where 

a similar problem was considered. The proof is completed in essentially 
the same way, i.e., via the definition of a suitable characteristic 


function. 
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APPENDIX IV 


P 
Reality of Aw of equation (5-47) 


It is assumed that T, and cf are real, as they are invariably simple 
known combinations of thermal averages. 


a 
From equation (5-51), 2. (é,) af are real forall 
P 


Let (R) = Real part of a 
Af (9) = Imaginary part of al 
4 
‘Then f f e 
= €,) Ax (A) +1 2G) A, (9) = reel »  (A-44) 
: r 
i.e. 


S EYAL (Neo bad (5 
po! 


Writing equation (A-45) in matrix form for { having the values 


Oy lites da gh@15 Fes 


t T* teva. t AL (9) 4 
a oo aes Az (4) o 
EVES .. ie) ee at 
! \n 

eG) A, (4) ‘ 


(A-46) 
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it is seen that the " es matrix", being a Vandermonde matrix and having 


all the Ep distinct, is non-singular, leading to-.the conclusion that 


A, (3) O 


A-47 
0 (A-47) 


all « . This completes the proof. 
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